Iterative Solution of Skew-Symmetric Linear Systems

نویسندگان

  • Chen Greif
  • James M. Varah
چکیده

We offer a systematic study of Krylov subspace methods for solving skew-symmetric linear systems. For the method of conjugate gradients we derive a backward stable block decomposition of skew-symmetric tridiagonal matrices and set search directions that satisfy a special relationship, which we call skew-A-conjugacy. Imposing Galerkin conditions, the resulting scheme is equivalent to the CGNE algorithm, but the derivation does not rely on the normal equations. We also discuss minimum residual algorithms, review recent related work, and show how the iterations are derived. The important question of preconditioning is then addressed. The preconditioned iterations we develop are based on preserving the skew-symmetry, and we introduce an incomplete 2× 2 block LDLT decomposition. A numerical example illustrates the convergence properties of the algorithms and the effectiveness of the preconditioning approach.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Convergence Properties of Hermitian and Skew Hermitian Splitting Methods

In this paper we consider the solutions of linear systems of saddle point problems‎. ‎By using the spectrum of a quadratic matrix polynomial‎, ‎we study the eigenvalues of the iterative matrix of the Hermitian and skew Hermitian splitting method‎.

متن کامل

A Preconditioner for Generalized Saddle Point Problems

In this paper we consider the solution of linear systems of saddle point type by preconditioned Krylov subspace methods. A preconditioning strategy based on the symmetric/ skew-symmetric splitting of the coefficient matrix is proposed, and some useful properties of the preconditioned matrix are established. The potential of this approach is illustrated by numerical experiments with matrices fro...

متن کامل

Numerical Investigation of Krylov Subspace Methods for Solving Non-symmetric Systems of Linear Equations with Dominant Skew-symmetric Part

Numerical investigation of BiCG and GMRES methods for solving non-symmetric linear equation systems with dominant skew-symmetric part has been presented. Numerical experiments were carried out for the linear system arising from a 5-point central difference approximation of the two dimensional convection-diffusion problem with different velocity coefficients and small parameter at the higher der...

متن کامل

Minimal Residual Methods for Complex Symmetric, Skew Symmetric, and Skew Hermitian Systems

While there is no lack of efficient Krylov subspace solvers for Hermitian systems, few exist for complex symmetric, skew symmetric, or skew Hermitian systems, which are increasingly important in modern applications including quantum dynamics, electromagnetics, and power systems. For a large, consistent, complex symmetric system, one may apply a non-Hermitian Krylov subspace method disregarding ...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:
  • SIAM J. Matrix Analysis Applications

دوره 31  شماره 

صفحات  -

تاریخ انتشار 2009